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A model of coupled maps with Pareto behavior 
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Abstract. A deterministic system of coupled maps is proposed as a model for economic activity among 
interacting agents. The values of the maps represent the wealth of the agents. The dynamics of the system 
is controlled by two parameters. One parameter expresses the growth capacity of the agents and the 
other describes the local environmental pressure. For some values of the parameters, the system exhibits 
nontrivial collective behavior, characterized by macroscopic periodic oscillations of the average wealth of 
the system, emerging out of local chaos. The probability distribution of wealth in the asymptotic regime 
shows a power law behavior for some ranges of parameters. 
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The study of wealth distribution in western societies 
has been a focus of much attention in the emerging area of 
research of Econophysics. A power law behavior is found 
in many economic activities, for instance, in income dis- 
tributions. Typically, high income earners, amounting to a 
few percent of the population, are distributed following a 
Pareto-like distribution. [1,2,3,4,5 . Hence, inequality in 
the wealth distribution is a fact in most economic activi- 
ties. The origin of such behavior seems to be caused by the 
interaction of the macro with the microeconomy. Here we 
propose a simple deterministic spatiotemporal model for 
economic dynamics where inequality emerges as a result 
of the dynamical processes taking place only at the micro- 
scopic scale. That is, the microeconomy fully determines 
the macroeconomic characteristics of the system. 

The model [B] consists of N interacting agents repre- 
senting companies, countries or other economic entities, 
placed as nodes on a network. For simplicity, we shall as- 
sume that the agents are distributed on a one-dimensional 
lattice with periodic boundary conditions. The state of an 
agent i, i — 1---N, is characterized by a real variable 
xt(i) £ [0, oo] denoting its wealth or richness at the dis- 
crete time t. The system evolves in time synchronously. 
Each agent updates its state xt{i) according to its present 
state and the states of its nearest neighbors. Thus, in a 
first approach, we propose that the value of Xt+i (i) is given 
by the product of two terms; the natural growth of agent 
i given by r(i)xt(i) with positive local ratio r(i), and a 
control term that limits this growth with respect to the 



local field 



**(*) = -[»*(» + + 1)], 



(1) 



through a negative exponential with parameter a(i), 

x t+ x[i) = r(i) x t {i) exp(- | x t (i) - a(i)& t (i) ). (2) 

The parameter r{i) represents the capacity of agent i to get 
richer and the parameter a(i) describes the local selection 
pressure [7] . This means that the largest rate of growth for 
agent i is obtained when x(i) ~ a(i)Wt(i), i.e., when the 
agent has reached some kind of adaptation to the local 
environment. In this paper we consider a homogeneous 
system with a uniform capacity r and a fixed selection 
pressure a for all the agents. 

The simplest collective behavior of the one-dimensional 
coupled map system described by Eqs. ([l} and ([2]) cor- 
responds to the synchronized or spatially uniform state. 
This state satisfies Xt(i) = x t and tf't(i) = x f , Vi, and its 
evolution is determined by the single map 



x t +i = rx t exp(- | (1 - a)x t ). 



(3) 
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For r < 1 the synchronized dynamics Eq. Q relaxes to 
the stable fixed point x = 0, while for r > 1 the synchro- 
nized system displays a sequence of bifurcations through 
different periodic and chaotic attractors, except for the 
singular case a = 1. For r > 1 the stable fixed point 
is x = logr/ 1 — a\. This point becomes unstable by a 
flip bifurcation at r = e 2 . For increasing r, the entire 
period-doubling cascade characteristic of unimodal maps 
and other complex dynamical behaviors are generated by 
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Eq. (|3|). However, it can be shown that such synchronized 
states are unstable. When a perturbation is introduced in 
the initial uniform state, the asymptotic dynamical state 
of the system is found to be more complex. 

In general, the collective behavior of the system can 
be characterized through the instantaneous mean field of 
the network, defined as 

H t =-J2Mj)- (4) 

3=1 

Figures [TJ a) andQJb) show a bifurcation diagram of H t as 
a function of the parameter r, for two different values of 
the selection pressure a. In Fig. QTa) it can be seen that H t 
reaches stationary values with some intrinsic fluctuations 
due to the local chaotic dynamics. However, for some pa- 
rameter values, nontrivial collective behavior [9] can arise 
in this system as shown in Fig.QJb). In this case, a macro- 
scopic variable such as the mean field or average wealth 
in the system follows a periodic behavior coexisting with 
chao ' " ' ' ' 1 
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Fig. 1. Ht as a function of r, for two values of a. For each 
value ofr, H t is plotted for 100 iterations, after discarding 9900 
transients. System size N = 10 4 . (a) a = 0.67. (b) a = 0.1. 

The statistical properties of the system can be ex- 
pressed in terms of the probability distribution of wealth 
among the agents. For some values of the parameters r 
and a, the distribution of wealth displays a power law be- 
havior, as shown in Fig. 2. For the values of parameters 
chosen, the probability distribution scales as P(x) ~ s a , 



with a = —2.86, similarly to scaling behaviors of Pareto 

type Hirprflir nlitninprl frnm flf+nfil prnnnrnv Hfttn [&] 
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Fig. 2. Log-log plot of the asymptotic probability distribution 
of wealth P(x) vs. agent wealth x. The distribution P(x) is 
calculated averaging the outcomes att — 10 4 of 100 realizations 
of random initial conditions. Parameter values are r = 12, 
a = 0.67, TV = 10 4 . The slope obtained is a = —2.86 with a 
correlation coefficient (5 — 0.98. 

In summary, we have shown that a power law behavior 
in the probability distribution of wealth can arise for some 
values of parameters in a deterministic system of interact- 
ing economic agents, such as in the coupled map model 
considered here. In addition, nontrivial collective behav- 
ior, where macroscopic order coexists with local disorder, 
can emerge in this system. 
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